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ABSTRACT 

The supernova remnant G347.3-0.5 emits a featureless power-law in X-rays, thought 
to indicate shock-acceleration of electrons to high energies. We here produce a broad- 
band spectrum of the bright NW limb of this source by combining radio observations 
from the Australia Telescope Compact Array (ATCA), X-ray observations from the 
Advanced Satellite for Cosmology and Astrophysics (ASCA), and TeV 7-ray observa- 
tions from the CANGAROO imaging Cerenkov telescope. We assume this emission is 
produced by an electron population generated by diffusive shock acceleration at the 
remnant forward shock. The nonlinear aspects of the particle acceleration force a con- 
nection between the widely different wavelength bands and between the electrons and 
the unseen ions, presumably accelerated simultaneously with the electrons. This allows 
us to infer the relativistic proton spectrum and estimate ambient parameters such as 
the supernova explosion energy, magnetic field, matter density in the emission region, 
and efficiency of the shock acceleration process. We find convincing evidence that the 
shock acceleration is efficient, placing > 25% of the shock kinetic energy flux into rel- 
ativistic ions. Despite this high efficiency, the maximum electron and proton energies, 
while depending somewhat on assumptions for the compression of the magnetic field in 
the shock, are well below the observed 'knee' at ~ 10 15 eV in the Galactic cosmic-ray 
spectrum. 
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1. INTRODUCTION 

The supernova remnant G347.3-0.5 (ROSAT catalog source name: RX J1713. 7-3946) is one 
of a growing number of supernova remnants (SNRs) showing nonthermal X-ray emission attributed 
to shock accelerated electrons. 5 In some cases, the nonthermal emission is evidenced by continuum 
emission extending to ~10 keV or above, usually superimposed with weak lines, or, as with SN1006, 
G347.3-0.5, and G266.2-1.2, by a complete lack of line emission in the observed X-ray band. 6 In all 
cases where the emission is not associated with a compact object, the nonthermal X-rays have been 
interpreted as synchrotron emission from shock accelerated TeV electrons (e.g., Mastichiadis & de 
Jager 1996; Allen et al. 1997; Tanimori et al. 1998; Reynolds & Keohane 1999). If this is true, the 
existence of these TeV electrons has important implications, not only for the production of cosmic 
rays, but also for the thermal properties of the shock heated, X-ray emitting gas in SNRs. 

Here we compare broad-band observations of G347.3-0.5 with results from a model combining 
hydrodynamic simulations of SNR evolution (e.g., Blondin & Ellison 2001, and references therein) 
with nonlinear diffusive particle acceleration (e.g., Ellison, Berezhko, &: Baring 2000). The observa- 
tions combine new radio data from the Australia Telescope Compact Array (ATCA) with existing 
X-ray (Slane et al. 1999) and 7-ray data (Muraishi et al. 2000). Our broad-band comparison of 
observations and model fits is a powerful way to constrain model parameters, but uncertainties 
from directly comparing radio, X-ray, and 7-ray observations are unavoidable and should be kept 
in mind. In an attempt to limit errors from uncertainties in emission volumes, we only consider 
observations from the northwest portion of G347.3— 0.5 where the center of the TeV 7-ray response 
function coincides with the brightest radio and X-ray emission. 

Previous work (e.g., Reynolds 1998) indicates that the steepness of the nonthermal X-ray 
synchrotron spectrum in shell-like SNRs results because these photons are generated by electrons 
in the steeply falling, high-energy tail of the distribution. The actual maximum electron energy, 
however, depends on the magnetic field strength, B2, in the emission region. If these electrons 
are accelerated by the first-order Fermi shock acceleration mechanism as we assume, it is almost 
certain that the same shocks accelerate ions, but no clear evidence of superthermal ions has yet 
been obtained (i.e., no pion-decay feature has been unambiguously identified at ~ 100 MeV in 
SNRs). Until a pion-decay feature is observed, all information on the energetic ions must be 
obtained through modeling. In shocks accelerating particles efficiently, the unseen ions dominate 
the shock dynamics and largely control the electron emission features (e.g., Drury 1983; Jones & 
Ellison 1991). 



5 As of this writing, the list of SNRs with nonthermal X-ray emission associated with shock acceleration rather 
than, or in addition to, pulsar-powered emission includes: Tycho, Kepler, Cas A, SN1006, G266.2-1.2, G347.3-0.5, 
RCW 86, G156.2+5.7, and 3C397. 

6 SN1006 shows line emission from the central regions, but nearly featureless spectra from the bright rim (e.g., 
Koyama et al. 1995). 



-3- 



Our main results are: (i) A spherically symmetric, wind-shell model with efficient acceleration 
can produce broad-band continuum spectra in excellent agreement with the observations of G347.3- 
0.5; (ii) Test-particle models of G347.3— 0.5 (i.e., ones where less than a few percent of the total 
energy ends up in relativistic particles) can be excluded with good confidence; (iii) The maximum 
electron energy in G347.3-0.5 is well below 10 14 eV. Thus, G347.3-0.5 joins several other SNRs 
now known not to produce electrons above 10 14 eV (e.g., Reynolds & Keohane 1999), well below 
the so-called knee in the Galactic cosmic-ray spectrum at ~ 10 15 eV. The maximum proton energy 
depends somewhat on whether or not the ambient magnetic field is compressed by the shock. If 
the field is compressed as the density, the maximum proton energy is approximately the same as 
the electron E max . In the other extreme, where no compression of the field occurs, we find that 
G347.3-0.5 could currently be producing protons to ~ 70 TeV and Fe +26 to ~ 10 15 eV. Unless 
substantial magnetic field amplification occurs, G347.3— 0.5 will not produce higher energy protons 
as it ages; and (iv) Our main conclusions concerning maximum particle energies and acceleration 
efficiency, are relatively insensitive to uncertainties stemming from the morphology of G347.3— 0.5. 

There are strong arguments for a SNR origin of cosmic rays based mainly on ion composition 
(e.g., Ellison, Drury, <fe Meyer 1997; Meyer, Drury, & Ellison 1997) and total energy requirements 
(e.g., Axford 1981; Blandford &: Eichler 1987; Drury, Markiewicz, & Volk 1989). The observations 
of nonthermal X-ray emission showing that TeV electrons are produced in SNRs adds support for 
this argument. However, there is still no direct evidence that ions are generated in particular SNRs 
in a fashion consistent with cosmic-ray observations. The low maximum proton energy inferred 
for G347.3-0.5 may present a problem for the scenario that SNRs are the main sources of Galactic 
cosmic rays below the knee and suggest that the cosmic-ray proton component above ~ 10 14 eV 
might come from some subset of SNRs substantially different from G347.3— 0.5. Alternatively, 
young SNRs which only produce ~ 10 13 eV particles when they are bright in X-rays, may somehow 
be able to produce higher energy particles at later times, perhaps by amplifying the magnetic field 
(e.g., Keohane, Gotthelf, & Petre 1998; Jun & Jones 1999; Lucek & Bell 2000). We note that the 
high magnetic fields inferred for Cas A make it likely that this particular remnant does produce 
~ 10 15 eV protons (e.g., Allen et al. 1997). 



2. Observations 

G347.3— 0.5 was observed with the Advanced Satellite for Cosmology and Astrophysics (ASCA) 
on 1997 March 25 for a total of 50 ks using three pointings of 10 ks duration and one (in the fainter 
southeast region) of 20 ks duration (Slane et al. 1999). Unlike the line-dominated X-ray spectra 
typical of SNRs, the spectrum for G347.3— 0.5 is observed to be featureless and to extend to at 
least 10 keV, where the ASCA effective area becomes very small. A power law spectrum provides 
an excellent fit to the X-ray data, with some indication that the spectral index varies between 
the brighter shell regions and the more diffuse interior and eastern regions. The brightest X-ray 
emission occurs along the northwest limb of the SNR (see Figure 1), where there is some indication 
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that the remnant may be interacting with a molecular cloud or a stellar wind bubble shell, although 
the low inferred density for this region more likely indicates some density enhancement caused by 
a perturbation of the pre-SN wind. The unabsorbed flux density from this region of the remnant 
is F x (0.5 — lOkeV) = 1.6 x 1CT 10 erg cm _2 s _1 with a power law photon index of 2.411°'^ . 
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Fig. 1.— G347.3-0.5 at 1.36 GHz, from the ATCA, with X-ray contours from the ROSAT PSPC 
(Pfeffermann & Aschenbach 1996). The modeling presented here corresponds to emission from the 
northwest limb of the SNR. The radio flux used for our modeling is from the two bright filaments 
indicated with white arrows. The X-ray spectrum was extracted from the region indicated by a box. 
The circle on the figure is placed at the centroid of the 7-ray emission detected by the CANGAROO 
telescope; the radius of the circle corresponds to the telescope point response function. 

Described here are new radio observations of G347.3-0.5 made with the Australia Telescope 
Compact Array (ATCA), a six-element synthesis telescope located near Narrabri, NSW, Australia. 
Observations at 1.36 GHz were carried out in January, March, and April 1998, in each case in 
a different array configuration to give maximal u — v coverage. The field was imaged in a 10- 
point mosaic, to ensure that the entire extent of the SNR, as well as nearby bright sources, were 
in the field-of-view. The interferometric visibilities were edited to remove corrupted data and 
interference. Calibrations were then applied to solve for the time-variable atmospheric gains above 
each antenna and for the polarization leakage in the receivers. An absolute flux density scale was 
applied to the data using observations of PKS B1934-638, for which a flux density of 15.0 Jy was 
assumed at a frequency of 1.4 GHz. The visibilities from all 10 pointings were then inverted to 
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produce a single image, which was deconvolved using the mosaiced maximum entropy approach 
described by Sault, Staveley-Smith, & Brouw (1996). The final image was then produced by 
applying a Gaussian restoring beam of dimensions 46" x 36". The sum of the flux from the two 
bright filaments in the northwest region of the remnant (indicated in Figure 1 and labeled as arc 1 
and 2 in Slane et al. 1999), corresponding to that from which the X-ray flux above was measured, 
is 5(1. 36 GHz) = 4 ± Uy. As noted in Slane et al. (1999), while arc 2 (the right-most arc in 
Figure 1) is askew from the perceived SNR shell, it does overlap an X-ray filament. Arcs 1 and 2 
have approximately the same radio flux. 

G347.3-0.5 has also been detected at TeV energies with the CANGAROO 3.8m atmospheric 
imaging Cerenkov telescope (Muraishi et al. 2000). As with the X-ray and radio observations, the 
~ 5.6a detection is from the NW rim of the remnant (indicated with a circle in Figure 1) and has 
an estimated flux of 5.3 ± 0.9[statistical] ± 1.6 [systematic] x 10 -12 photons cm~ 2 s _1 at energies 
> 1.8 ± 0.9 TeV. 

The fluxes from these three detections are shown in Figures 3 and 4 where we have increased 
the observational error bars on the radio and TeV points by a factor of two (light-weight bars) to 
compensate for uncertainties in the relative emission volumes. In particular, the increased radio 
error bar will accommodate the possibility that arc 2 is not associated with the SNR shell and thus 
none of our results will change if this is the case. 

3. Model 

3.1. Nonlinear diffusive shock acceleration 

Collisionless shocks are known to accelerate particles efficiently. In the heliosphere, where 
shocks are observed directly with spacecraft, there is clear evidence that the quasi-parallel earth 
bow shock (with a sonic Mach number, M$o < 10) can place 10-30% of the solar wind kinetic energy 
flux into superthermal particles (e.g., Ellison, Mobius, & Paschmann 1990). Interplanetary shocks 
(IPSs) also accelerate ambient particles but generally with lower efficiencies due to their lower Mach 
numbers (generally Mgo < 3 for IPSs) (e.g., Baring et al. 1997). However, on at least two occasions, 
exceptionally strong IPSs have been observed to accelerate particles with high efficiency (Eichler 
1981; Terasawa et al. 1999). Hybrid plasma simulations showing direct injection and acceleration 
of ambient particles are consistent with these observations (e.g., Scholer, Trattner, & Kucharek 
1992), as are convection-diffusion models (e.g., Kang & Jones 1997). 

Compared to shocks in the heliosphere, it might be presumed that the strong forward and 
reverse shocks in young SNRs are much more efficient. As far as we can tell, there are no important 
differences in typical plasma parameters (i.e., density, magnetic field, temperature, composition) 
between the solar wind and the ambient ISM. The only real differences come from the much larger 
shock spatial and time scales and much higher flow speeds and Mach numbers in SNRs. In general, 
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these factors should increase the acceleration efficiency substantially over heliospheric shocks, but 
some qualifications must be made. 

At quasi-parallel shocks (those where the upstream magnetic field makes an angle, @Bn, of less 
than about 45° with the shock normal), injection from the shocked heated plasma should occur eas- 
ily and the fraction of incoming ram kinetic energy flux going into superthermal particles increases 
with Mach number and shock size and/or age (as measured in particle gyroradii or gyroperiods). 
If however, the fraction of shocked heated particles injected is small enough, high Mach number, 
test-particle solutions can result (see Berezhko Sz Ellison 1999). In quasi-perpendicular shocks (ones 
where 45 < 0B n ^ 90°), injection and acceleration of shock heated thermal particles is suppressed 
above some Mach number unless strong turbulence is present (see Ellison, Baring, & Jones 1996, 
for example). Highly oblique interplanetary shocks are observed to accelerate the ambient solar 
wind (even without pickup ions, Ellison, Jones, & Baring 1999), but these are low Mach number 
shocks. As the Mach number increases, injection becomes more likely if the turbulence is strong 
enough for cross-field diffusion in the downstream plasma to be important. 

While the shock obliquity will surely vary around a SNR and some portions of the shocks 
may be quasi-perpendicular with low injection rates, other regions should be parallel enough for 
efficient injection and acceleration to occur. Radio emission is a sure sign that relativistic electrons 
are present and flux levels, sharp radio edges, etc., show that these particles are locally accelerated in 
young SNRs rather than just compressed Galactic cosmic-ray electrons (e.g., Achterberg, Blandford, 
& Reynolds 1994). The shock acceleration and evolutionary models we use here ignore effects from 
oblique magnetic fields, other than those that are mimicked by a low injection efficiency. We do 
assume, however, that the magnetic field responsible for synchrotron radiation can be compressed 
by the shock. For a recent discussion of how the injection efficiency may vary with 0B n in SNRs, 
see V61k (2001). 

The nonlinear effects in efficient shock acceleration are of three basic kinds: (i) the self- 
generation of magnetic turbulence by accelerated particles, (ii) the modification (i.e., smoothing) 
of the shock structure by the backpressure of accelerated particles, and (iii) the increase in shock 
compression ratio as relativistic particles are produced and some high-energy particles escape from 
the system. Briefly, (i) occurs when counter-streaming accelerated particles produce turbulence 
in the upstream magnetic field which amplifies as it is convected through the shock. This am- 
plified turbulence results in stronger scattering of the particles, and hence to more acceleration, 
quickly leading to saturated magnetic turbulence levels near 5B/B ~ 1 in strong shocks. This is 
the so-called Bohm limit where the mean free path is of the order of the gyroradius, and in this 
limit, the distinction between parallel and oblique shocks blurs. The wave-particle interactions 
produce heating in the shock precursor which may be observable and which lowers the overall ac- 
celeration efficiency. Effect (ii) occurs because the accelerated population presses on the upstream 
plasma and slows it. An upstream precursor forms, in which the flow speed (in the shock frame 
and in the absence of instabilities) is monotonically decreasing. Since particle diffusion lengths 
are increasing functions of momentum, high momentum particles sample a broader portion of the 
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flow velocity profile, and hence experience larger effective compression ratios than low momentum 
particles. Consequently, higher momentum particles have a flatter power-law index than those at 
lower momenta (producing a concave upward spectral curvature) and can dominate the pressure 
in a non-linear fashion. We call the ratio of densities spanning the entire upstream precursor the 
overall compression ratio, r tot , to distinguish it from the small scale subshock compression ratio, 
r su b, where most of the shock heating occurs. For a more complete discussion of the nonlinear 
effects of diffusive shock acceleration, see Berezhko & Ellison (1999). 

The conservation of mass, momentum, and energy fluxes, i.e., the Rankine-Hugoniot relations, 
yield the standard expression for the compression ratio in shocks where the magnetic field is not 
dynamically important and an insignificant fraction of energy is placed in superthermal particles 
(i.e., test-particle shocks): 

_ ( 7 + l)Mg 
rtot " ( 7 - 1)M| + 2 ' W 



where M^o = y Po ^jk/(7-Fb) is the sonic Mach number, 7 is the ratio of specific heats, po is the 
unshocked density, Po the unshocked pressure, and is the shock speed. 7 In test-particle shocks, 
rtot = r suo . The ratio of shocked to unshocked pressure follows: 

P 2 _ 2 7 M| Q - (7 - 1) 

P ~ 7 + 1 ' { ' 

If M so » 1, 

rtot ^ ^ , (3) 
7-I 

and 

F2 = Fo 2^ = ^. (4) 
7+1 7+I 

For a electron-proton plasma where the downstream temperatures are equal, the post-shock proton 
temperature in this test-particle limit, T tp , is: 

tP - r to t fc( 7 + 1) ' { ) 

where k is Boltzmann's constant and m p is the proton mass. 

Young SNR shocks are expected to have high Mach numbers and if the test-particle situation 
is assumed with 7 ~ 5/3, we expect r to t — 4. However, when shocks accelerate particles the overall 
compression ratio increases above test-particle values (even in non-radiative shocks) for two reasons 
(effect iii above). First, as relativistic particles are produced and contribute significantly to the total 
pressure, their softer equation of state makes the shocked plasma more compressible (7 — > 4/3). 
Second, and most important, is that particles escape from strong shocks draining energy flux which 



Here and elsewhere, the subscript '0' indicates unshocked values and the subscript '2' indicates shocked values. 
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must be compensated for by ramping up the overall compression ratio to conserve the fluxes. Just 
as in radiative shocks, this is equivalent to 7 — > 1. For strong injection, it can be shown (Berezhko 
& Ellison 1999) that r to t is an ever increasing function of Mach number, i.e., 



where Mao is the Alfven Mach number (see also, Kazanas & Ellison 1986; Malkov 1997). Simul- 
taneously with rtot increasing, the subshock compression ratio, r su b, which is mainly responsible 
for heating the gas, must drop below the test-particle value. Thus, high compression ratios are ac- 
companied by low post-shock temperatures (in efficient shocks, the post-shock proton temperature 
can easily be 1/10 T tp ). 9 If injection is weak enough, however, equation (6) does not apply and 
high Mach number shocks do not accelerate enough particles to cause the total compression ratio 
to increase above ~ 4 (see Berezhko & Ellison 1999, for the parameters describing this condition). 

The most important parameters associated with nonlinear shock acceleration are the Mach 
numbers (i.e., the shock speed, V s ^, pre-shock hydrogen number density, and preshock magnetic 
field, Bo), the injection efficiency, 7?inj,p (i-e., the fraction of total protons which end up with su- 
perthermal energies), and the maximum proton energy produced, E max>p . As described in Berezhko 
& Ellison (1999), our model includes Alfven heating in the precursor which reduces the efficiency 
compared to adiabatic heating. 

We refer the reader to Ellison, Berezhko, k, Baring (2000) for a discussion of our treatment of 
electron injection and acceleration. 



Slane et al. (1999) suggest that the forward shock of G347.3— 0.5 is still within the pre-SN 
stellar wind bubble with part of it interacting with denser material. We describe this evolution 
with a spherically symmetric, wind-shell model with a pre-SN density distribution as shown in the 
top panel of Figure 2. The initial density profile in the ejecta has a power law distribution, p oc r~ n , 
with a constant density plateau at small radii, as described in Chevalier (1982a). Since we only 
consider the forward shock at fairly late stages in the SNR evolution, our results are insensitive to 
n and the position of the plateau. They are also insensitive to whether a power-law or exponential 
ejecta mass density distribution (e.g., Dwarkadas & Chevalier 1998) is assumed. For concreteness, 



8 As long as the speeds of the escaping particles are much greater than the shock speed, escaping momentum and 
mass fluxes can be neglected (Ellison 1985). 

9 Hughes, Rakowski, & Decourchelle (2000) have used this effect to reconcile the high values of the shock speed 
determined from expansion measurements and the low post-shock electron temperature determined from X-ray line 
models in SNR E0102. 2-7219. 




if 1 < Mf < M A0 ; 
if 1«M A0 <M S 2 , 



(6) 



3.2. 



SNR Evolution and Model Parameters 
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Fig. 2. — The top panel shows a snapshot of a typical density profile at the start of the hydrody- 
namic simulation (i.e., at t snr = 0). The spatial scale assumes that G347.3— 0.5 has a diameter of 
~ 45' and is 6 kpc away. The inner edge of the shell is at ~ 34 pc and the current radius of the 
forward shock is at ~ 40 pc. The middle panel shows the evolution of the forward shock speed as a 
function of radius and the bottom panel shows the shock speed as a function of remnant age. The 
different curves in the lower two panels were obtained by running the hydro code with the indicated 
effective ratios of specific heats. In all cases, the speed drops sharply as the forward shock enters 
the shell. 
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we use n = 9 in all of our models. As in Chevalier (1982b), the density of the pre-SN wind is given 
by 

p w = A r~ 2 , (7) 

where Aq = dM\ oss /dt/(4:TTv w ), dM\ oss /dt is the pre-SN mass loss rate, and v w is the constant 
wind speed. At some distance from the explosion, we place a shell of uniform density, p s heii = 
m p ra s heii (1+4 f a ) — our spherically symmetric representation of the density enhancement suggested 
by Slane et al. (1999). Here, f a = 0.1 is the number fraction of helium and n s h c ii is the proton 
number density in the shell. While the radio image of G347.3— 0.5 suggests a shell with the strongest 
emission coming from the NW rim, the complex morphology makes an exact determination of the 
diameter somewhat ambiguous. Some radio and X-ray emission extends across ~ 60', but a smaller 
~ 30' diameter ring-like structure is also present. As a compromise, we adopt an angular diameter 
of 45', and to further limit parameters, we perform all of our fits assuming a distance of D snr = 6 
kpc (Slane et al. 1999). With these assumptions, the outer shock is at a radius i? s k ~ 40 pc and 
we take (somewhat arbitrarily) the inner radius of the shell to be -Ri n ~ 34 pc, giving a shell 
thickness approximately equal to the thickness of radio arc 1, i.e., (-R s k — i?i n )/i? s k ~ 3'/22.5'. 
Fortunately, none of our conclusions depend critically on these particular values. Furthermore, as 
long as p w (R- m ) "C p s hcii, our results are insensitive to Aq and we use Aq = 6xl0 12 gm cm -1 which 
is consistent with a wind extending to 34 pc. 10 

The evolution is calculated with the Virginia Hydrodynamics (VH-1) time-dependent hydro- 
dynamic numerical code in one dimension where it is implicitly assumed that dynamic effects from 
magnetic fields are negligible. For a description of this code, see Blondin & Ellison (2001) and 
references therein. With the shock speed and other parameters obtained with the hydro code, we 
calculate the electron, proton, and helium spectra and with them the synchrotron, bremsstrahlung, 
inverse-Compton, and pion-decay photon continuum spectra. We then compare the photon emission 
to the observations and vary parameters until a satisfactory fit is obtained. 



3.2.1. Approximations of the wind-shell model 

Of necessity, our model is approximate in some important ways. First, we assume the SNR 
is spherically symmetric when, in fact, the emission is strongest from one section of the rim where 
the shock is likely interacting with denser material. We compensate for this at the most basic level 
by defining the emission volume, V^, m i s , as a fraction, / vo i, of the shell volume encompassed by the 
forward shock, i.e., 




10 For dM loBB /dt = 7xl(T 6 M Q yr -1 and v w = 60 km s _1 (giving A = 6xl0 12 gm cm -1 ), ~ 4M Q would exist in a 
spherical wind extending to 34 pc. With this v w , the wind bubble would have taken ~ 6x 10 5 yr to form and could 
still not have stagnated in a low density interstellar medium. 
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From Figure 1 we estimate that 0.05 < / vo i < 0.25. 

Second, we use a one-zone model where the particle spectra and photon emission are assumed 
to be generated by a shock of constant speed during the time the forward shock is in the shell. 
We ignore any evolution of the spectra or the shock during the time the shock is in the shell and 
the shock speed used for acceleration is the average speed in the shell as determined by the hydro 
code. Furthermore, we neglect any contributions from the reverse shock or reflected shocks which 
occur when the forward shock enters the shell. These should be excellent approximations since 
the forward shock speed doesn't vary greatly while in the shell (middle panel Figure 2), particles 
produced earlier in the low density wind or ejecta have been reduced by adiabatic losses, and 
reflected shocks in the shell are weak compared to the forward shock. 

The evolution of the SNR will be influenced by the acceleration process through the increased 
compression ratios of the forward and reverse shocks and by the escape of high energy cosmic rays. 
We have approximated this effect by performing our hydrodynamic calculations with values of the 
effective ratio of specific heats, 7 e ff < 5/3, as described in detail in Blondin & Ellison (2001). The 
effect on the forward shock speed is shown for typical parameters in the bottom panel of Figure 2. 
The solid curve is for j c g = 5/3 (r to t — 4), the dashed curve is for 7 e g- = 1.4 (r to t — 6), and 
the dotted curve is for 7 e g = 1.1 (r to t — 20). When 7 c g- = 5/3, the shock reaches 40 pc after 
t snr ~ 1.3 x 10 4 yr, for 7 c g- = 1.4, it reaches 40 pc after ~ 1.45 x 10 4 yr, and for j e g = 1.1, 40 pc is 
reached after ~ 1.8 x 10 4 yr. While globally changing 7 e g only approximates the effects of particle 
acceleration on the evolution (particle escape is not explicitly included and the evolution of 7 e g is 
not modeled, for example), we do not believe errors from this approximation are significant since 
all of the models we discuss below have compression ratios less than ~ 6 when the forward shock 
is in the shell. As indicated in Figure 2, at -R s k = 40 pc V s \r varies by less than 25% for r to t < 6. 
Despite this small change in V^k, the particle spectra produced by a shock with r to t = 6 are still 
strongly nonlinear and differs importantly from test-particle spectra as shown in the fits below. 

Finally, our model ignores X-ray line emission which, for a thermal plasma, can be much 
more intense than the pure bremsstrahlung component. For example, the 0.2-10 keV flux for a 
kT = 0.2 keV Raymond-Smith model (Raymond & Smith 1977) is ~ 20 times higher than for a 
thermal bremsstrahlung model with the same temperature and emission measure. This discrepancy 
decreases with temperature, dropping to ~ 4 at kT = 1 keV and ~ 1.4 at kT = 2.5 keV. If the 
featureless X-ray spectrum in G347.3— 0.5 is synchrotron emission, the underlying bremsstrahlung 
continuum which we calculate must be low enough so lines, had they been calculated, would not 
appear. 

3.2.2. Model parameters 

Most of the parameters used in our models are described in Ellison, Berezhko, & Baring (2000) 
and/or shown here in Table 1. Briefly: (i) The explosion energy, E sn , largely determines the overall 
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dynamics and is typically assumed to be near 10 erg. The ejecta mass, M e j, is unimportant for 
the ages and parameters we consider as long as it is less than about 3 Mq. If M c j > 3M Q , the 
reflected shock produced when the forward shock enters the shell can interact with the density peak 
behind the reverse shock producing complicated flow structures. To avoid these complications, we 
limit consideration to M c j < 3M and do not list this parameter in the Table, (ii) The magnetic 
field affects the particle acceleration and largely determines the shape of the broad-band photon 
emission through its influence on synchrotron emission. The field morphology is certain to be 
complex in any real SNR and will vary with location (see Jun & Jones 1999, for an excellent 2- 
D MHD study of shock-cloud interactions). To consider a range of possibilities in our idealized, 
spherically symmetric model, we assume the downstream field, B2, is either equal to the upstream 
(unshocked) field, i.e., B2 = Bq, or that B2 is compressed along with the gas, i.e., B2 = r tot Bq. The 
actual situation should lie somewhere in between. We assume the unshocked proton temperature, 
T p Q = 10 4 K = T e Q = T a Q, where T e o (T a o) is the upstream electron (helium) temperature. Our 
solutions are insensitive to the upstream temperature as long as it is below ~ 10 6 K. The shocked 
electron to proton temperature ratio, T C 2/T P 2, is a free parameter which has not yet been deduced 
from first principles for collisionless shocks. The electron temperature is an important parameter 
for thermal bremsstrahlung and X-ray line emission, but is not strongly constrained for G347.3— 0.5 
because no lines are seen. In any case, we assume T e 2/T p 2 = 1 in all of our models but note that if 
electron-proton equilibration timescales are long compared to dynamic timescales, T e 2/T p 2 could be 
considerably less than one, lowering the thermal bremsstrahlung to synchrotron ratio. 11 (iii) The 
injection efficiency, f]mj,p, is the fraction of thermal protons that end up as accelerated particles. 
This parameter influences the overall acceleration efficiency, (iv) The electron to proton ratio at 
relativistic energies, (e/p) re i, is an arbitrary parameter in our models but is expected to lie between 
0.01 and 0.05 if (e/p) re i in Galactic cosmic rays is typical of that produced by the strong shocks in 
young SNRs. This factor influences the normalization of the photon emission from electrons and is 
important for determining the pion-decay contribution to TeV 7-rays relative to inverse-Compton. 
The maximum energy cosmic rays obtain depends on the scattering mean free path, A, which is 
assumed to be, 

= Vmip r g,max (j"g/ r g,max) j (9) 

where ry m f p is a constant, r g = p/ (qB) is the gyroradius in SI units, r giIliax is the gyroradius at 
the maximum momentum, p max , and a is a constant parameter. Small values of r/ m f p imply strong 
scattering and allow higher maximum proton energies in a given system. We use the Bohm limit 
(^mfp = 1) m a U °f our models and note that while good fits to the observations can be obtained 
with larger r? m f p , these models will have similar acceleration efficiencies with lower proton maximum 
energies than the models we show, (v) Finally, n s h e ii is the proton number density in the shell which 
models the density enhancement presumed to exist in the NW rim. Slane et al. (1999) conclude, 



The combination of high post-shock density and low post-shock proton temperature means that electrons should 
equilibrate on a faster timescale in shocks undergoing efficient acceleration compared to those where little acceleration 
occurs (Decourchelle, Ellison, & Ballet 2000). 
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from the lack of thermal emission, that the mean density around the remnant is quite low. With 
large uncertainties, they estimate upper limits on the ambient density of 0.014-0.28 (D snr /6kpc) -1 / 2 
cm -3 , with a higher upper limit (< 1 (D sm /6 kpc) -1 / 2 cm~ 3 ) for the northwest rim. 

The output values are also given in Table 1. The shock speed, V s k, is the average obtained from 
the hydro simulation over the shell thickness between R = 34 and 40 pc. The age of G347.3— 0.5 is 
unknown, but Slane et al. (1999) have used emission measures to estimate a range 2000 < t snr < 
40, 000 yr. The value in the Table is that necessary to yield an outer shock radius of i? s k ~ 40 
pc consistent with an angular diameter ~ 45' at a distance of 6 kpc. The time the shock spends 
between 34 pc and 40 pc is At^. Given V^k and the density, temperature, and magnetic field in 
the shell, the sonic and Alfven Mach numbers are determined. The above parameters allow the 
nonlinear shock model to determine r tot , r su b, and 7 e ff, as well as the proton, electron, and helium 
spectra. Note that consistent values of r to t and j e s are found by iterating the hydro simulation 
with the particle acceleration calculation. The shocked magnetic field is E>2, and the shocked proton 
temperature is T p 2. The shocked proton temperature for the corresponding test-particle shock, T tp , 
is also given. 

The particle momentum distributions, f(p), are calculated as in Ellison, Berezhko, Sz Baring 
(2000) with the turnover at the highest energies described by: 



f(p) -> f(p) exp 



1 ( P 



a \p 



(10) 



Here, a is the same as in Eqn. (9) and p max is determined by setting the acceleration time equal to 
At s h, or by setting the diffusion length of the highest energy particles equal to 1/10 R s k, whichever 
gives the lowest p max (see Baring et al. 1999). It turns out that a relatively slow turnover (i.e., 
a ~ 0.5) is required to obtain the proper X-ray slope and to join the radio and X-ray observations. 
A similar effect was required to match the intensity and shape of the X-ray observations in SN1006 
(Reynolds 1998; Berezhko, Ksenofontov, & Petukhov 1999; Ellison, Berezhko, & Baring 2000), and 
we take a = 0.5 in all of our models. The electron and proton energies corresponding to p max 
are -E m ax,c and -E m ax,p respectively and are listed in the Table. The acceleration efficiency, e re i, is 
defined as, 

£rcl " (l/2) P0 Vl ' (Uj 

where F re i and F esc are the fractions of energy flux in relativistic protons and escaping protons, 
respectively. 

The spectral index at 2 GHz, where F u cx v a , is shown in the next to last row of Table 1. Radio 
data has been collected at 2.4 GHz from G347.3— 0.5 but has not yet been analyzed. When the 
2.4 GHz flux is combined with the existing 1.4 GHz flux, the inferred spectral index can be used 
as a check of our models and may help constrain parameters. Finally, the last row in the Table is 
the fraction of the shell volume, / vo i, producing the emission seen from the NW limb. Based on 
equation (8), models with / vo i > 0.25 are unacceptable. 
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Table 1: Forward shock parameters for G347.3— 0.5. All models listed have r/ m f p = 1, a = 0.5, 
T p0 = 10 4 K, and L> snr = 6 kpc. 



Input 


A 


B 


C 


E SQ [10 51 erg] 


2 


3 


4 


Bo \mG] 


2.4 


2.4 


16 


Vinj,p 


lxlO" 3 


5xl0~ 6 


1 x 10~ 3 


(e/p)re\ 


0.03 


0.05 


0.04 


n s he\\ [cm -3 ] 


0.015 


0.02 


0.008 


Output 


V sk [km s- 1 ] 


980 


1100 


2100 


tsm [yr] 


9000 


8000 


5000 


At sh [yr] 


5000 


4700 


2200 


M so 


60 


70 


130 


M A0 


27 


35 


6.3 


not 


6.1 


4.1 


4.3 


^sub 


3.0 


4.0 


3.4 


7cff 


1.4 


1.65 


1.6 


B 2 \mG] 


15 


9.6 


16 


kT p2 [keV] 


0.33 


1.2 


3.4 


kT tp [keV] 


1.0 


1.2 


4.5 


^max,p [TeV] 


14 


16 


70 


£max,c [TeV] 


14 


16 


32 




0.55 


0.03 


0.20 


ct(2GHz) 


-0.45 


-0.51 


-0.54 


Flux 


/vol 


0.16 


6.8 


0.19 
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-12 -9 -6 -3 3 6 9 
Log 1Q E y [MeV] 

Fig. 3. — Comparison of broad-band emission from G347.3-0.5 with diffusive shock acceleration 
models. Model A uses a likely set of parameters with B2 = r tot Bq, model B is a fit obtained with a 
test-particle shock, and model C was obtained with parameters chosen to obtain the highest proton 
energy. In all panels, the short-dashed curves are synchrotron, the dot-dashed are bremsstrahlung, 
the dotted are inverse-Compton, the long-dashed curves are pion-decay, and the normalization is 
set by matching the synchrotron fluxes at ~ 5 keV. The heavy line at X-ray energies is the fit to 
the NW rim observations from Slane et al. (1999). The gamma-ray point is from Muraishi et al. 
(2000), and the radio data is first published here. The light-weight error bars on the radio and TeV 
points are included to compensate for the possible difference in relative emission volumes between 
the various energy bands. In the top panel we show the expected continuum sensitivity limits for 
the GLAST spacecraft (Gehrels k Michelson 1999), the INTEGRAL spacecraft (Winkler 1996), 
and the VERITAS array (Weekes et al. 1999). The solid curve in B is the total emission with 
/vol = 0.25. 
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3.3. Fits to G347.3-0.5 

The broad-band emission from G347.3— 0.5 allows us to constrain the parameters of the shock 
model to fairly narrow ranges. In the top panel of Figure 3 (model A), we show a satisfactory fit 
to the broad-band observations with B2 = r tot BQ. The parameters for this model and two others 
discussed below are listed in Table 1. The four components of the emission are labeled and it's clear 
that the synchrotron emission produces the radio and X-ray emission, while the TeV point is fit 
with inverse-Compton photons from electrons scattering off the 2.73 K background radiation. The 
value / V oi = 0.16 is consistent with the morphology and n s h e ii = 0.015 cm -3 is consistent with the 
upper limit estimate of Slane et al. (1999). We note that using a lower E sn would force a smaller 
n s heii and a larger / vo i; E sn > 1 x 10 51 erg is clearly favored. 

Using model A as a reference, we investigate how sensitive the fit is to particular parameters. 
In Figure 4 we show variations in Bq, n s h e ii, E sn , and r/i n j,p- In each panel we keep all parameters 
except the one labeled fixed to the values used in model A. For all models, we have normalized the 
synchrotron emission to the flux in the X-ray band and plot the sum of the four photon components. 
The solid curve in each panel is model A. In the top panel, the dotted curve is for Bq = 1 /iG, the 
solid curve is for Bq = 2.4 (iG, and the dot-dashed curve is for Bq = 7 fj,G. In all cases, the number 
in parentheses is f vo \. The dependence on Bq is quite strong and the shape of the broad-band 
emission clearly sets upper and lower limits on Bq. Furthermore, the Bq = 1 /xG result can be 
excluded because f vo \ ~ 4; a low B produces too little flux. 

The changes produced by varying n s h e ii are shown in the second panel of Figure 4. The dotted 
curve is for n s h c ii = 0.005 cm~ 3 , the solid curve is for ra s heii = 0.015 cm~ 3 , and the dot-dashed curve 
is for n s hcii = 0.04 cm -3 . Clearly, the n s h e ii = 0.005 cm~ 3 result is too flat in the X-ray band, while 
the n s heii = 0.04 cm~ 3 result is somewhat too steep. The value f vo \ ~ 1.1 for n s h e ii = 0.04 cm~ 3 , 
shows that this high density example produces too little flux. This occurs because is low in this 
case (520 km s _1 ) producing a low -E max ,e and forcing a higher normalization to match the high 
energy X-rays. 

In the third panel, we vary E sn from 0.5 x 10 51 erg (dotted curve) to 2x 10 51 erg (solid curve) 
to 4xl0 51 erg (dot-dashed curve). The low explosion energy gives an unacceptable fit because of 
the shape in the X-ray band and because / vo i ~ 8. However, for E sn = 4 x 10 51 erg the shape is 
still reasonably good considering the uncertainties in the model and the normalization (/ vo i ~ 0.03) 
while low, cannot be excluded. 

In the bottom panel of Figure 4 we show results varying the injection efficiency. The dotted 
curve is for r^i n j iP = 10 -5 , the dot-dashed curve is for ?7i n j,p = 10 -4 , and the solid curve is model 
A with r?inj,p = 10~ 3 . The shape of the sum of the emission is virtually identical for ?7inj,p = 10~ 4 
and 10~ 3 . However, the intensity of the emission varies with rj- m j tP and for r/i n j,p = 10~ 4 , / vo i ~ 0.3, 
while for ??mj,p = 10~ 3 , / vo i ~ 0.16. The constraint that the emission is only observed from a small 
fraction of the shell implies that ??inj,p > 10 -4 . This will be true unless E sn is considerably greater 
than 2 x 10 51 erg. The ??mj,p = 10~ 5 result, which is what would be expected if the shock did not 
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-12 -9 -6 -3 3 6 9 



Log 10 E y [MeV] 

Fig. 4. — These plots and data are as in Figure 3 except here the curves are the sum of synchrotron, 
bremsstrahlung, inverse-Compton, and pion-decay emission. In each panel, the labeled parameter 
is varied as all other parameters are held fixed with the values shown in Table 1 for model A. The 
numbers in parentheses are f vo \. These models all have B2 = r to t -Bo- I n the bottom panel, / vo i ~ 10 
for ??mj,p = 10~ 5 (dotted curve), f vo \ ~ 0.3 the ??mj,p = 10 -4 (dot-dashed curve), and / vo i ~ 0.16 for 
Vinj,p = 10 -3 (solid curve). 
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accelerate particles efficiently, is unacceptable because of the overall shape and low emissivity (i.e., 

/vol - 10). 

To further investigate the limits on the acceleration efficiency, we choose r/inj.p low enough so 
only a few percent of the total energy flux ends up in relativistic particles and freely adjust other 
parameters to obtain a fit. With ??mj,p = 5 x 10 -6 (yielding e rc i ~ 3%) we are able to obtain an 
approximate match to the shape of the observed emission with Model B in Figure 3 and Table 1. 
Note that this model requires E sn = 3 x 10 51 erg; lower energies give poorer fits. While the shape 
of the continuum emission is a reasonable match to the observations, the normalization requires 
/vol ~ 6.8, well above an acceptable value. Furthermore, the synchrotron emission from the test- 
particle shock is somewhat too steep to simultaneously match the radio and X-ray fluxes, and the 
ratio of synchrotron to bremsstrahlung continuum emission in the X-ray range is less than a factor 
of 3. If we further decrease r?i n j,p 5 the synchrotron to bremsstrahlung ratio will drop and at some 
point X-ray emission lines would stand above the continuum flux in conflict with the observations. 
This conclusion depends, of course, on T e 2/T p 2 and could be weakened by assuming T e 2/T p 2 <C 1. 
However, the relative acceleration efficiency of electrons and protons is also likely to scale as T^/T^ 
(Ellison, Berezhko, & Baring 2000) making it difficult to produce the observed radio and TeV flux 
if T C 2/T p 2 is too small. Furthermore, if a flatter synchrotron slope is required to give a better match 
between the radio and X-rays (as in model ^4), this forces the shock to be more efficient (to produce 
a larger rtot) and implies a larger r/mj )P . The overall flux, the lack of X-ray emission lines, and the 
radio to X-ray flux ratio all point toward efficient shock acceleration. Using this model we can 
get a firm lower limit on the acceleration efficiency by increasing ??mj,p until / vo i ~ 0.25. We find 
/vol i$ 0.25 for r/inj,p > 3x 10~ 5 , yielding a lower limit e re i > 0.25, i.e., even with E sn = 3x 10 51 erg, 
the shock must place at least 25% of its kinetic energy flux in relativistic particles. Note that we 
took (e/p) re i = 0.05 for this example. A lower (e/p) TC \ would imply a higher acceleration efficiency. 

In order to investigate the maximum proton energies that G347.3— 0.5 can generate, we have 
produced Model C (bottom panel of Figure 3) where we have chosen parameters (including taking 
B2 = Bo) with the intention of maximizing -E max , p . Very generally, diffusive shock acceleration 
yields the same maximum energy per charge for all particles (e.g., Drury 1983; Baring et al. 1999), 
so the only way the proton and electron spectra can have different maximum energies is for losses 
(synchrotron and/or inverse-Compton) to be important for the electrons. Unless we have neglected 
an abundant supply of photons in addition to the 2.73 K background radiation, inverse-Compton 
losses are not important and only when the magnetic field is large can -E max , P be much greater 
than £ maX)e . However, a large field produces a lower Alfven Mach number which weakens the 
acceleration and, in order to fit the broad-band spectrum, a higher shock speed is required. This 
can be achieved with a larger E sn and/or a lower n s h e ii, and we have taken E sn = 4x 10 51 erg 
and n s heii = 0.008 cm~ 3 for Model C. The observed radio flux sets the strength of the magnetic 
field in the emission region combined with the relativistic electron density, i.e., the flux is loosely 
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proportional to B 2 n e . Furthermore, the TeV inverse-Compton emission comes from the same 
electrons that produce the X-rays. Since the inverse-Compton emission is independent of B, it 
constrains the number of relativistic electrons. The combination of these constraints allows us to 
limit -Emax.p and we find in our most optimistic case that E maXyP < 70 TeV. 

In Figure 5 we show the momentum distribution functions, f(p), for models A, B, and C, 
where we have plotted the dimensionless quantity, p 4 f(p)/{m p cn p o), to emphasize the spectral 
curvature. In this representation, the test-particle result of f(p) oc p~ 4 is a horizontal line (model 
B). The concave shape of f(p) (at superthermal energies and below the cutoff) is evident for models 
A and C, as is the large difference between the relativistic particle populations in these shocks with 
efficient acceleration compared to the inefficient shock B. As noted above, a flat electron spectrum 
(at energies greater than m p c 2 ) produced by a large r to t ( r tot — 6 in model A) allows a good match 
between the radio and X-ray observations. 

The arrows in Figure 5 indicate the positions of -E m ax,p and E max>e as determined from equa- 
tion (10) with a = 0.5. In all cases, f(p) starts to fall off exponentially at momenta less than 
Pmax and much below p = 10 6 m p c. Furthermore, if G347.3-0.5 continues to evolve uniformly with 
no magnetic field amplification, it will not produce higher energy particles as it ages (for a more 
complete discussion of the time particles take to obtain their maximum energies, see Berezhko & 
Volk 1997, 2000). Since the causes of magnetic field amplification are still uncertain, it remains 
possible that higher fields and therefore higher energy particles could be produced in the future. 
Magnetic field amplification apparently occurs in Cas A (e.g., Gull 1975; Keohane, Gotthelf, & 
Petre 1998) but if it is due to Ray leigh- Taylor instabilities in this young remnant as commonly 
assumed, it is unlikely that this mechanism would occur in G347.3— 0.5, which is much older. 



4. Conclusions 

G347.3-0.5 is one of a growing number of shell-type SNRs showing nonthermal X-rays believed 
produced by synchrotron emission from TeV electrons accelerated by the remnant shocks. In an 
attempt to have a consistent, broad-band spectrum, we have combined recent X-ray and TeV 7-ray 
observations of the remnant with new radio observations and have restricted all of the observations 
to the bright northwest rim of the shell. Comparing a spherically symmetric, wind-shell model of 
SNR evolution and efficient diffusive shock acceleration against the observations shows that the 
observed emission can be well fit across ~ 18 decades in energy with a single set of parameters all 
close to expected supernova and ISM values. The densities we find for the emission region (~ 0.01 
protons cm -3 ) are consistent with the upper limits found earlier by Slane et al. (1999), and are 
below those expected for either a molecular cloud or a swept-up shell produced as the pre-SN wind 



This is not a strict proportionality because changing either B or the density produces changes in the acceleration 
model and/or the SNR evolution. 
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Fig. 5. — Proton (solid curves) and electron spectra (dashed curves) for models A, B, and C used 
to generate the photon spectra shown in Figure 3. The arrows indicate E maXtP and -Emax.e- The 
dotted curve in the bottom panel shows the high momentum portion of a possible Fe +26 spectrum 
from Model C. The normalization relative to protons is that shown in Figure 8 of Ellison, Drury, 
& Meyer (1997). 
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interacts with the interstellar medium (or with a slower red supergiant wind). This suggests the 
shock is still within the low density bubble but has encountered a density enhancement perhaps 
caused by perturbations in the pre-SN wind. 

We put limits on the acceleration efficiency and find that G347.3— 0.5 is currently putting at 
least 25% (and more likely ~ 50%) of the forward shock ram energy flux into relativistic particles 
(i.e., cosmic rays). Recently, Hughes, Rakowski, & Decourchelle (2000) have inferred efficient shock 
acceleration in SNR E0102. 2-7219 in the Small Magellanic Cloud from the shock speed, determined 
from expansion measurements, and the postshock electron temperature, determined from X-ray line 
models (unlike G347.3— 0.5, lines are visible above the synchrotron continuum in the post-forward 
shock gas). They concluded that the extremely low postshock temperature can only be reconciled 
with the high shock speed if a sizable fraction (~ 50%) of the shock energy goes into cosmic 
rays. Our results are fully consistent with this and with the important point that efficient particle 
acceleration causes a strong coupling between the broad-band emission from relativistic particles 
and the properties of the shocked heated X-ray emitting gas (see Decourchelle, Ellison, & Ballet 
2000). We have verified that Sedov models (not shown) with no density enhancement produce 
satisfactory nonlinear fits to the observed spectral shape with input parameters similar to those 
given here, but are inconsistent with the observed remnant diameter and distance. 13 We conclude 
that the nonlinear aspects we infer, in particular the required high acceleration efficiency, do not 
depend significantly on the geometric details of the wind-shell model while the overall normalization 
does. 

Furthermore, we find that the maximum energies of the cosmic-ray electrons and protons 
are well below 10 15 eV. While it has been known for some time that SNRs emitting nonthermal 
X-rays do not produce electrons to energies near 10 15 eV (e.g., Reynolds & Keohane 1999), the 
nonlinear effects we model make it possible to infer the proton distribution from photons emitted 
by electrons. For no magnetic field compression in the shock and other parameters adjusted to 
produce the highest possible E max ^ p and still give a reasonable fit to the observations, we find a 
proton maximum energy < 70 TeV. More realistic parameters with magnetic field compression yield 
E max , P ~ 20 TeV. However, SNR shocks will also accelerate heavier ion species such as He, C, O, 
and Fe, and the maximum energy in diffusive shock acceleration scales as charge. Therefore, this 
SNR could produce a power-law Fe +26 spectrum to ~ 10 15 eV. To illustrate, we have estimated 
the Fe +26 spectrum produced by our model C and placed it in the bottom panel of Figure 5. The 
normalization is determined by the cosmic abundance of iron (~ 3.1 x 10~ 5 iron nuclei for every 
proton) plus the enhancement expected because iron is injected and initially accelerated when it 
is trapped in dust grains. In fact, the normalization of iron relative to hydrogen in Figure 5 is 



' In standard Sedov models, the relatively high shock speeds and ambient densities which produce acceptable fits 
to the observed shape occur at shock radii which are too small to be consistent with a remnant diameter of 45' (or 
even 30') at D S m ~ 6 kpc. Only if D S nr < 3 kpc could standard models produce acceptable fits. The wind-shell 
model allows the shock to propagate for large distances in the low density wind before encountering the higher density 
emission region. 
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just that shown in Figure 8 of Ellison, Drury, & Meyer (1997) at 10 8 keV/nucleon. Cosmic-ray 
observations in the 10 14 — 10 16 eV range typically measure the total energy of the particles and have 
large uncertainties in differentiating individual ion species (e.g., Bernlohr et al. 1998; Glasmacher 
et al. 1999; Swordy & Kieda 2000). Nevertheless, a number of authors have produced models of 
this transition energy range adding together the fluxes from different ion species with cutoffs that 
scale as charge (e.g., Stanev, Biermann, & Gaisser 1993; Erlykin & Wolfendale 2000). Our high B$ 
results are in general agreement with this work. 
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